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ABSTRACT 

We present model lightcurves which have been created in order to explain the orbital 
modulation observed in the radio emission of Cyg X-l. We invoke variable absorption 
by the stellar wind as the black hole jet orbits around the OB companion star and 
find that a very simple model is able to reproduce the amplitudes and frequency 
dependence of the observed lightcurves. 

Key words: stars:individual:Cyg X-l — stars:winds, outflows — radio contin- 
uum:stars — binaries: general 



1 INTRODUCTION 

Cygnus X-l was first observed to be a radio source by Braes 
& Miley (1971). At the time of discovery the radio flux was 

, ~ 20 mjy and it generally maintains a constant mean level 
of ~ 10-25 mjy at cm wavelengths while the system is in 
its usual low/hard X-ray spectral state. However, when Cyg 

■ X-l makes one of its rare state changes to the soft state the 
radio emission appears to be quenched and drops to below 
a detectable level (Pooley 2001; Brocksopp et al. 1999a, and 
references within). This behaviour is similar to that of GX 
339—4 - another black hole X-ray binary which has been 
monitored in the radio throughout a soft state period and 
transitions (Fender et al. 1999, Corbel et al. 2000). 

An important feature of the radio emission is the ob- 
served flat synchrotron spectrum to millimetre wavelengths 
(Fender et al. 2000) with no evidence for a cut-off at higher 
frequencies. Unfortunately the infrared emission is domi- 
nated by the supergiant and so it is not possible to determine 
how far the flat spectrum continues. This spectrum suggests 
that Cyg X-l has a quasi-continuous jet of the type modelled 
by Hjellming & Johnston (1988) and/or Blandford & Konigl 
(1979). These predictions have recently been confirmed with 
the VLBI imaging of the jets of Cyg X-l (Stirling et al. 
2001). 

Two types of periodic behaviour can be observed in the 
radio emission. The longer of these is at ^140 days (although 
is probably not stable in the long term) and is possibly due 
to precession of the disc-jet system (Brocksopp et al. 1999a); 
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alternatively a variable mass accretion rate such as that pro- 
posed for LMC X-3 (Wilms et al. 2001; Brocksopp, Groot 
& Wilms 2001) is also a potential cause of this long period. 
The 5.6-day orbital period can also be detected in the radio 
emission (Pooley, Fender & Brocksopp 1999) and it is this 
modulation that we now investigate. 

The orbital period of Cyg X-l has been well-established 
for more than two decades and particularly since the work 
of Gies & Bolton (1982). The ephemeris has been revised 
recently (e.g. Brocksopp et al. 1999b, Sowers et al. 1998, 
LaSala et al. 1998) but the results have all been consistent 
with each other. 

X-ray data have also revealed the orbital period (Pa- 
ciesas et al. 1997) and in Brocksopp et al. (1999a) we dis- 
cussed the significance of finding the orbital period at all 
wavelengths in the Cyg X-l system. Detection of the or- 
bital period of black hole X-ray binaries in optical and in- 
frared photometry is standard and double-peaked orbital 
lightcurves, due to the gravitational pull on the star by the 
black hole, are expected - minima occur in the lightcurves at 
the two conjunctions. However, it is not so common to detect 
the orbital period in X-ray or radio data, particularly since 
Cyg X-l is a non-eclipsing system and is thought to have 
a circular orbit. Thus Brocksopp et al. (1999a) postulated 
that these modulations could be the result of absorption by 
the stellar wind. This was later confirmed by Wen et al. 
(1999) who studied the orbital modulation of the RXTE 
soft X-ray data and produced a relatively successful model 
of wind absorption. The X-ray orbital lightcurves have also 
been studied by Balucinska-Church et al. (2000) who found 
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that the duration and phase of the minima are also variable 
and likely to be due to 'dumpiness' in the stellar wind. 

At radio wavelengths we found that the modulation 
is frequency-dependent with the strongest modulation at 
higher frequencies - contrary to what one would expect for 
straightforward free-free absorption of a point source. There 
is also a phase lag present, the duration of which may also be 
frequency-dependent. These observations add support to the 
proposed continuous jet model on account of an extended 
jet-like structure probing different densities of an absorbing 
medium. Furthermore it would be unlikely that the modula- 
tions could be maintained as consistently as observed, were 
the radio emission produced in discrete ejections. 

We now produce a model of wind absorption that can 
explain the observed radio emission at three wavelengths. 
To do this we use the focussed stellar wind model of Friend 
& Castor (1982), which they and Gies & Bolton (1986b) 
have shown to be applicable to the Cyg X-l system. In this 
model the spherical stellar wind of the supergiant becomes 
distorted by the gravitational field of the black hole, thus 
resulting in enhanced emission occuring close to the axis 
joining the supergiant and black hole. 



2 OBSERVATIONS 

We have 15 GHz radio data from the Ryle Telescope of the 
Mullard Radio Astronomy Observatory, Cambridge, UK; de- 
tails of the observations can be found in Brocksopp et al. 
(1999a, and references therein). Further data have been ob- 
tained from the Green Bank Interferometer which monitored 
a number of sources at 2.25 and 8.3 GHz up until October 
2000 - details are given in Waltman et al. (1994). 

The data points at each frequency were taken between 
1996 November and 1999 June and have had the flaring 
episodes removed in order to smooth the lightcurves. They 
were then folded on the orbital period according to the 
ephemeris of Brocksopp et al. (1999b). 



3 MODEL LIGHTCURVES 

In Pooley, Fender & Brocksopp (1999) it was discovered that 
the radio emission of Cyg X-l was modulated on the or- 
bital period with an amplitude that, perhaps surprisingly, 
increased with the frequency of the emission; not what we 
would expect for a typical free-free absorption law. This in- 
dicates firstly that the density of the wind cannot be treated 
as constant along the line of sight - indeed, there is no stellar 
wind model which considers the wind density to be constant 
- and secondly that the optical depth along the jet is not 
constant. 

We use the work of Leitherer, Chapman & Koribalski 
(1995, and references within) to determine the (frequency 
dependent) radius from the supergiant at which the stellar 
wind becomes optically thin to free-free absorption of ra- 
dio emission (r ~ 0.244, measured from infinity; Wright & 
Barlow 1975): 

' * ' ' " (1) 



W v = 4.0 x W L 'T 



-1/2 1/3 



Rr. 



where we assume a wind temperature T ~ 5000K (assum- 
ing the results of e.g. White & Becker 1983 and Leitherer, 



Chapman & Koribalski 1995), use M ~ 2.0 x 10" 6 M /yr 
(Friend & Castor 1982) and ~ 1580 km/s (Gies & Bolton 
1986b) and calculate the Gaunt factor, g v (see Leitherer, 
Chapman & Koribalski 1995). 

The resultant values of W v are 4.6 x 10 14 , 1.8 x 10 14 
and 1.2 x 10 14 cm for v = 2.25, 8.3 and 15 GHz respectively. 
These values are considerably greater than the binary sepa- 
ration and so we must consider two different absorption com- 
ponents: phase-dependent absorption takes place within the 
binary orbit, providing the modulation, and additional (non 
phase-dependent) absorption takes place in the foreground 
between the orbit and the wind radius (which lowers the flux 
by a constant amount). 

We assume that the radio emission is absorbed by the 
wind via the free-free mechanism (Rybicki & Lightmann 
1979): 

- 4 (2) 



ft = 0.018T 3//2 z 2 n^riiU 2 gjf = 
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Then the observed flux is: 



(3) 



where So is the intrinsic flux density of the radio source 
and the absorption components, k$ (phase-dependent) and 
ft/ (foreground) are determined below. 

Unfortunately there is no way of knowing what the in- 
trinsic flux density of the radio source is - therefore we have 
used the observed flux density of 14.0 mjy ± 4.0, 6.0, 1.2 
mjy for the 2.25, 8.3 and 15 GHz observations respectively 
(i.e. a flat spectrum); this value is consistent with the peak 
flux density at each frequency, allowing for the errors of the 
GBI and Ryle telescopes. The value of the intrinsic flux den- 
sity is somewhat arbitrary but its variation leads to minimal 
change in the amplitude of the modulation and so merely has 
the effect of shifting the model lightcurve up or down. 

In order to estimate the distance along the jet at which 
each frequency becomes optically thin (r ~ 1), we note that 
the orbital modulation at 2.25 GHz is negligible. Thus we 
assume that the 2.25 GHz emission site is approximately co- 
incident with the wind radius calculated previously. We then 
use the relationship r„ oc v~ l , as determined by Blandford 
& Konigl (1979) for a compact flat spectrum jet, to esti- 
mate the emission sites of the 8.3 and 15 GHz flux. Thus 
the emission site is located at: 



r u = 4.6 x 10 1 



2.25 
///GHz 



(4) 



We note that this frequency dependence on distance 
along the jet has already been observed, the best example 
of which is GRS 1915+105 (Mirabel et al. 1998). 



3.1 Phase-dependent absorption 

The focussed stellar wind model of Gies & Bolton (1986b) 
and Friend & Castor (1982) is used. Since we are interested 
in the radio emission (which originates a long way from the 
orbital plane) we use the parameters derived for the region 
> 20° from the binary axis. We take Equation 2 of Gies & 
Bolton (1986b) for the wind density: 



P<t> 



R 



Po 



(l-r*/RY 



(5) 
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Figure 1. Schematic showing the Cyg X-l system as considered in this simple model for the radio modulation. Note that the jets extend 
considerably further than shown in this diagram. The 'edge' of the wind is defined as where the optical depth of the wind (measured 
from infinity) drops to r ~ 0.244 (Wright & Barlow 1975). 



where r* is the radius of the supergiant (~ 1.387 x 10 12 
cm) and a is a constant (=1.05) as defined in Gies & Bolton 
(1986b). This is substituted into the equation for k above: 

6.4 x W 60 R- 4 u- 2 
^ ~ (1 - 1.387 x lO 12 ^- 1 ) 21 ( ' 

Thus the phase-dependent absorption can be determined 
at any point R(x, y, z) where x — y = z = corresponds to 
the centre of the OB star. 



The foreground absorption was actually found to be 
negligible (r/ < 10~ 36 between the 'edge' of the wind and 
the binary orbit) and all significant absorption took place 
as the black hole passed behind the stellar companion - this 
is in agreement with Gies & Bolton (1986a) who determine 
that the density of the wind beyond the binary orbit is very 
small. We note that recombination of the ionised wind at 
large distances from the supergiant has not been taken into 
account - therefore we might have expected to twer-estimate 

Kf. 



3.2 Foreground absorption 

Similarly we determine the foreground absorption for any 
value of R. However, once outside the orbit the equation 
of Gies & Bolton (1986b) is no longer valid - there it is 
sufficient just to consider the stellar wind as a spherical shell 
expanding at the value of Voo as far as W v . 



Pf 



Sttv 



3A V 2 



Again, this is substituted into Eq. 2 for k above: 



K f = 2.83 x W 1J u~ 2 (R' 



(7) 



(8) 



3.3 Results 

With knowledge of the absorption at each point in the stel- 
lar wind and the position of the radio emission site, it was 
possible to integrate along the line of sight from the observer 
to the circle traced by the jet about the star - i.e. a circle 
centred on x = 0, y — 0, z — r„. This was repeated for 360 
positions around the circle in the frame of the binary sys- 
tem. It was assumed that the angle between the line of sight 
and the circle traced by the jet remained constant and was 
equivalent to the inclination angle of the binary system. It 
was also assumed that all radio emission comes from a single 
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Figure 2. The radio data at 2.25, 8.3 and 15 GHz folded on the 
orbital period . Theoretical curves are also plotted showing that 
the simple model fits the data reasonably well. 



Table 1. Values of the reduced x 2 between the data and the 
model at each of the four inclination angles plotted in Fig. j^. The 
model fits reasonably well, particularly at 15 GHz for which the 
data quality is significantly better. 
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Figure 3. Predicted amplitudes for the orbital modulation at a 
range of radio frequencies. The model suggests that the emission 
is saturated at superior conjunction of the black hole for high 
frequencies and we hope to test this observationally with SCUBA 
observations. 



jet - Doppler boosting of the approaching jet and Doppler 
de-boosting of the receding jet lead to significant differences 
between the observed flux of each and so this is indeed a 
reasonable assumption. 

The resultant model lightcurves are shown in Fig. ^| 
plotted over the folded data. Each plot shows four differ- 
ent curves, corresponding to a range of inclination angles. 
We note that while the curves appear sinusoidal the peaks 
are in fact slightly broader than a sinusoid - this is in 
agreement with previous attempts to fit a sine wave to the 
mean lightcurves which were somewhat unsuccessful (Poo- 
ley, Fender & Brocksopp 1999) . The reduced \ 2 between the 
data and the model has been determined and the results 
listed in Table |. 

Since the model has been defined to fit the 2.25 GHz 
data with zero modulation there is not much that can be 
said about the top plot in Fig. ^. However, its lack of orbital 
variability shows that we were justified in considering the 
2.25 GHz emission to be produced outside the wind radius. 
The 8.3 and 15.0 GHz lightcurves are reasonably well- fit by 
the model, particularly at an inclination angle of 40° which 
is consistent with values quoted in the literature. Clearly it 



would be beneficial to improve on the sensitivity of the 2.25 
and 8.3 GHz data in order to confirm the quality of the fit. 
We would, however, like to add a note of caution regarding 
the use of our work to predict the inclination angle of the 
system - the focussed stellar wind model of Gies & Bolton 
(1986b), on which our model is largely based, assumes an 
inclination ~ 30° and so to suggest that our model predicts 
i ~ 40° is not necessarily appropriate without further devel- 
opment of the wind model as suggested below. 

We can also use our model to predict the amplitude 
of the orbital modulation at radio frequencies other than 
those observed here. The results are shown in Fig. ^| for 
four values of the inclination angle and it is clear that the 
degree of absorption at superior conjunction of the black 
hole increases with frequency. Once the frequency has been 
increased to ~50 GHz the model predicts saturation at an 
orbital phase of ~ 0.15 and this may be expected given that 
the high frequency emission is produced right at the base of 
the jet in the most dense regions of the stellar wind - this 
is yet to be tested observationally and we have proposed 
millimetre observations with which to do this. 
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4 DISCUSSION 

4.1 Improving the model 

The most obvious problem with this basic model is that it 
has been assumed that the radio emission at each frequency 
is essentially a point source at position r„ along the jet, and 
that this position lies within the stellar wind. Since this as- 
sumption results in reasonable theoretical curves it is tempt- 
ing to consider this to be the case. However the VLBA maps 
of Stirling et al. (2001) show that, while it is a reasonable 
assumption for the 15 GHz emission, it is unlikely to be the 
case at 8.3 GHz - at each epoch there is a significant fraction 
of the jet well outside the wind radius and it is clearly not 
a point source. Re-running the code to allow for ~ 50% of 
the emission to be outside the modulated region produces 
too little modulation compared with the GBI data-set. Cer- 
tainly it would be beneficial to obtain higher sensitivity data 
in order to check this discrepancy. Alternatively it is possi- 
ble that the wind density in the outer parts of the wind has 
been underestimated and that the foreground absorption is 
actually significant - a higher intrinsic flux density would 
allow for < 100% of the emission to be modulated and to 
still produce sufficiently deep minima. Problems with the 
wind model are discussed below. 

There are a number of further assumptions which have 
been made in order to create the model lightcurves - we 
have not taken account of the wind structure, other than the 
simple assumption that we can use the focussed wind model 
as described above. Firstly it is unlikely that the wind is a 
constant temperature throughout - with the stellar surface 
thought to be at a temperature of ~ 30000 K there will 
be a region of hotter and less dense wind near the star. 
This may not affect the radio emission (at least at lower 
frequencies) since it comes from a region a long way from the 
star. Secondly we have assumed a constant density function 
throughout the wind. Again this is unlikely as the wind will 
be clumpy and irregular. Adding these corrections to the 
model at this stage is well beyond the scope of this paper 
- detailed modelling of the stellar atmosphere and wind in 
the Cyg X-l system is planned but is a long-term project. 

A further effect on the stellar wind is caused by the pres- 
ence of the X-ray source. The stellar wind of an OB star, 
such as the optical counterpart to Cyg X-l, is emitted by 
radiation pressure in the star's outer atmosphere which scat- 
ters and absorbs photons in various spectral lines formed in 
the wind. The acceleration of the wind is influenced by the 
degree of X-ray heating and ionisation from a compact ob- 
ject situated in the path of the wind (Hatchett & McCray 
1977). While an increase in X-ray flux increases the line- 
driving force and hence the wind velocity, once the X-ray 
luminosity is of the order 10 34 ergs/s the line force and wind 
velocity begin to decrease again. This is due to the increasing 
proportion of highly ionised species (emitting at optical, ul- 
traviolet and X-ray wavelengths), which leads to a decrease 
in luminosity of certain emission lines and subsequent al- 
teration of the wind spectrum^ (Hatchett, Buff & McCray 

t An effect of this appears to be observed — when the X-ray source 
of Cyg X-l makes a transition to the high/soft state (i.e. an in- 
creased source of soft X-ray photons at the centre of the accretion 
disc), a quenching of the Ha focussed wind emission is detected 



1976). By the time the X-ray luminosity reaches 5 x 10 34 
ergs/s the high ionisation in the region around the compact 
object is sufficient to overcome the line-driving force, thus 
quenching the radiative acceleration in a region known as 
the Stromgren zone (Macgregor & Vitello 1982). With a lu- 
minosity of ~ 10 37 erg/s, the X-ray source of Cyg X-l should 
produce a significant effect on the wind. However, assuming 
that the X-ray flux decreases as ~ 1/r 2 and that the X-ray 
ionisation is only significant in the high velocity wind close 
to the orbital axis (Hatchett, Buff & McCray 1976), it is 
probable that this Hatchett-McCray effect is insignificant at 
the radius of the radio emission sites. (See van Loon, Kaper 
& Hammerschlag-Hensberge 2001 for detailed modelling and 
size estimates of the Stromgren zone.) 



4.2 Further considerations 

An additional property that is not represented by the model 
is the phase offset seen in all radio lightcurves - the theo- 
retical curves have been shifted in phase by 0.15 in order 
to match the observations. While Pooley, Fender & Brock- 
sopp (1999) and Brocksopp et al. (1999b) find a frequency 
dependence to this offset, it is not confirmed in the ex- 
tended datasets used here. (The reason for the disppearance 
of the frequency-dependence to this lag is uncertain - it is 
of course possible that it was some artefact caused by the 
high noise levels in the GBI data. Higher sensitivity VLA 
observations are planned in order to investigate this effect 
further.) However, the offset is still very much present in the 
mean lightcurves. A possible cause of this could be some sort 
of drag effect - as the black hole passes through the stellar 
wind the jets may be 'swept back' and trail behind slightly. 
Frequencies emitted towards the far end of the jet (i.e. lower 
frequencies) would then be delayed by a greater amount than 
frequencies emitted nearer the black hole. 

Such a phase offset can be simply modelled - assuming 
that the phase delay is dominated by self-absorption of the 
jet, the jet will lag the black hole by a distance 

r v 

X Vorb 

Vjet 

where Vj e t and v or b are the jet and orbital velocities respec- 
tively. This corresponds to a phase lag of 

r v 

Vjet X Porb 

where P or b is the orbital period. By assuming ri5 G H Z ~ 6.8 x 
10 13 cm then a phase lag of 0.15 could be produced by a jet 
travelling at Vj e t ~ 0.03c - while this is mildly relativistic, 
it is probably an underestimate, particularly given the lower 
limits estimated by Stirling et al. (2001). 

As mentioned previously, the presence of the orbital 
period in the radio emission is unexpected - indeed the 
only other systems for which this has been the case are 
LSI +61°303 and Cir X-l (Taylor & Gregory 1982; Haynes, 
Lerche & Murdin 1980). In both cases the period is much 



via equivalent width measurements (Voloshina, Lyuty & Tarasov 
1997). We note that the radio jet is associated with the low/hard 
state only (Brocksopp et al. 1999a) and so this observation is not 
relevant to the model presented here. 



© ?? RAS, MNRAS 000, ??-§ 



6 Brocksopp, Fender, Pooley 



longer than that of Cyg X-l (26.5 and 16.6 days respec- 
tively) and the nature of the radio emission is also very dif- 
ferent. Rather than the smooth sinusoid of Cyg X-l, these 
two sources show periodic radio outbursts thought to be pro- 
duced at periastron of an eccentric orbit. Therefore even if 
the companion stars emit strong stellar winds (as certainly 
is the case for LSI +61° 303), it is unlikely that the wind 
absorption model can be applied here - it would predict a 
minimum at periastron as this is the point at which the wind 
is most dense. Alternatively the wind model may still be rel- 
evant but the effects are dominated by the higher accretion 
rate at periastron and subsequently greater jet power. How- 
ever, it is interesting to note that the only previous attempt 
to model the radio modulation of Cyg X-l is Han (1993). She 
also invokes an elliptical orbit with an inclined jet - this can 
be discounted given the lack of evidence for anything other 
than a circular orbit in optical spectroscopy presented in the 
literature (e.g. Brocksopp et al. 1999b). 

One important application and test of this absorption 
model will be the lightcurve prediction of jet emission at 
other frequencies, in particular the millimetre and infrared 
emission close to the black hole. Cyg X-l has been detected 
at millimetre wavelengths but there is insufficient data avail- 
able to detect any orbital modulation. A comparison of the 
model and lightcurves at higher frequencies would also be 
useful in probing the regions of higher ionisation around the 
black hole. 



5 CONCLUSIONS 

We have presented a model to explain the orbital modulation 
observed in the radio emission of Cyg X-l. By invoking ab- 
sorption of the radio emission by the stellar wind of the com- 
panion star we are able to produce model lightcurves with 
amplitudes and frequency dependence comparable with the 
real data. We acknowledge that the fits could be improved 
with the use of a more sophisticated stellar wind model - 
the production of such a wind model is work in progress. 
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